Localized Energetics of Linear Gravity: Theoretical Development 
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We recently developed a local description of the energy, momentum and angular momentum car- 
ried by the linearized gravitational field, wherein the gravitational energy-momentum tensor displays 
positive energy-density and causal energy-flux, and the gravitational spin-tensor describes purely 
spatial spin [1, 2]. We now investigate the role these tensors play in a broader theoretical context, 
demonstrating for the first time that (a) they do indeed constitute Noether currents associated 
with the symmetry of the linearized gravitational field under translation and rotation, and (b) they 
are themselves a source of gravity, analogous to the energy-momentum and spin of matter. To 
prove (a) we construct a Lagrangian for linearized gravity (a covariantized Fierz-Pauli Lagrangian 
for a massless spin-2 field) and show that our tensors can be obtained from this Lagrangian using 
a standard variational technique for calculating Noether currents. This approach generates for- 
mulae that uniquely generalize our gravitational energy-momentum tensor and spin tensor beyond 
harmonic gauge: we show that no other generalization can be obtained from a covariantized Fierz- 
Pauli Lagrangian without introducing second derivatives in the energy-momentum tensor. We then 
construct the Belinfante energy-momentum tensor associated with our framework (combining spin 
and energy-momentum into a single object) and as our first demonstration of (b) we establish that 
this Belinfante tensor appears as the second-order contribution to a perturbative expansion of the 
Einstein field equations, generating the gravitational field in a manner equivalent to the (Belinfante) 
energy-momentum tensor of matter. By considering a perturbative expansion of the Einstein-Cartan 
field equations, we then demonstrate that (b) can be realized without forming the Belinfante ten- 
sor: our energy-momentum tensor and spin tensor appear as the quadratic terms in separate field 
equations, generating gravity as distinct entities. Finally, we examine the role of field redefinitions 
within these perturbative expansions; in contrast to our tensors, the Landau-Lifshitz tensor is found 
to require a non-local field redefinition in order to be cast as a source of the gravitational field. In 
an appendix, we also give a brief treatment of the global quantities that our framework defines, and 
verify their equivalence (within the quadratic approximation) to the ADM energy-momentum and 
angular momentum. 



I. INTRODUCTION 

There are at least three ways to quantify the energy 
of a physical system. One approach is to consider in- 
teractions with a second system (the energy of which is 
known) and to seek a function-of-state of the first sys- 
tem which, by undergoing equal and opposite changes, 
accounts for the energy lost or gained by the second sys- 
tem. Alternatively, a Lagrangian for the physical system 
might be constructed, and the energy identified as the 
Noether charge associated with the Lagrangian's symme- 
try under translation in time. Thirdly, and most simply 
of all, one can "weigh" the system; the energy is then 
determined by the gravity it generates. 

In the preceding articles [1, 2], we arrived at a local 
description of the energy, momentum, and angular mo- 
mentum of the linearized gravitational field. 1 The re- 
sulting gravitational energy-momentum tensor T ab and 
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To clarify: we have not performed the demonstrably impossible 
feat of finding an energy-momentum tensor r a j, ~ V/iV/i and 
spin tensor s a bc ~ hV/i that are invariant under the linearized 
gauge transformations of the gravitational field Sh a i, = V( a £;,). 
Rather, we rely on a gauge-fixing programme (motivated by key 



spin tensor s a bc are particularly notable in that, when- 
ever the field is transverse-traceless, they describe non- 
negative energy-density, causal energy flow, and spatial 
spin; moreover, these properties, and the gauge invariant 
energetics of an infinitesimal probe, motivate a natural 
gauge-fixing procedure. These two tensors were derived 
by what is essentially the first method described above: 
we sought functions of the gravitational field which could 
account for the energy, momentum, and angular momen- 
tum exchanged locally with matter. The main purpose of 
this present investigation is to demonstrate that T a b and 
s a bc fulfil the other two roles of energy-momentum and 
spin: (a) they are Noether currents of local translations 
and rotations, and (b) they appear as self-interaction 
terms in the field equations, generating gravity alongside 
the energy-momentum and spin of matter. 

The paper is organized as follows. In section II, we 
construct a Lagrangian for linear gravity (a covariantiza- 
tion of the Fierz-Pauli Lagrangian for a massless spin-2 
field) and show that it generates T ab and s a bc according 
to standard variational definitions of energy-momentum 



properties of T ab and 



and the energetics of an infinitesimal 



gravitational detector) to remove the freedom to perform such 
transformations, and hence arrive at a physically unambiguous 
description. 
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and spin. This process of "gauging" the translational 
and rotational symmetries of the background, and de- 
riving T a b and s a bc as functional derivatives of the La- 
grangian with respect to the gauge fields, confirms their 
status as Noether currents of translational and rota- 
tional symmetry. (More precisely, the vectors r^ a are 
the Noether currents of translations in the x M direction, 
and j^J 1 = 2x[ ll T„] a + s a flL/ are the Noether currents of ro- 
tations in the x^a^-plane.) In section III we then demon- 
strate that, under a local redefinition of the gravitational 
field, and s a bc appear (combined into a single Be- 
linfante energy-momentum tensor) as the quadratic part 
of the vacuum Einstein field equations. The same tech- 
niques are then applied to Einstein-Cartan gravity, with 
Tab and s a bc appearing as the quadratic parts of two sep- 
arate field equations. In both cases, T a b and s a bc quantify 
the self-interaction of the gravitational field, and gener- 
ate gravity in an identical fashion to material energy- 
momentum and spin. Finally, in section IV we exam- 
ine the significance of the field redefinitions used in sec- 
tion III; in contrast to our tensors, we show that a non- 
local redefinition is required in order to cast the Landau- 
Lifshitz tensor [3] as a source of the gravitational field. 

From all of this we learn that T a b and s a bc stand on an 
equal footing with the energy-momentum and intrinsic 
spin of matter: they can be derived from the symmetries 
of a suitable Lagrangian, and behave as sources for the 
gravitational field. These developments solidify the ten- 
sors' physical interpretation, and embed them within the 
same theoretical apparatus that has been used to define 
gravitational energy- momentum in the past [3-6]. Fur- 
thermore, through their role in the non-linear field equa- 
tions, we gain insight into how T a b and s a bc may be ex- 
tended beyond the linear regime, and also uncover a new 
and possibly valuable set of gravitational field- variables. 

Our notation and conventions have not changed since 
the previous papers [1, 2]. We will, however, introduce 
a new variety of index: Greek letters with overbars, 
which we will use to enumerate the components of ten- 
sors in the non-holonomic basis defined by the tetrad 
e?; sec appendix A for details. (These bars should not 
to be confused with those placed over 2-index tensors, 
which indicate trace- reversal: h a b = h a b ~ fjabh/2 and 
Tab = Tab — QatiT /2.) As before, when working in the flat 
background spacetime (M,g a b) it will often be conve- 
nient to express our tensors in a coordinate system {x^} 
that is Lorentzian with respect to the background metric: 
9p.v = "H^v- 



II. LAGRANGIAN FORMULATION 



from a Lagrangian L that also generates the field equa- 
tions of linear gravity. 2 Before we do this, however, we 
will have to make some accommodation for the harmonic 
gauge condition, 

= 0, (2) 

which we have been diligently enforcing since its ap- 
pearance as an unexpected consequence of the deriva- 
tion of T^ [1]. While this condition has been immensely 
valuable, naturally reducing the gauge ambiguity of our 
framework, it now has the potential to interfere with the 
arbitrary field variations that occur when taking func- 
tional derivatives of a Lagrangian. To avoid this problem, 
we shall temporarily relax the gauge condition (2) and 
aim to derive from C an energy-momentum tensor and 
spin tensor that generalize r M „ and s a beyond the har- 
monic gauge, reducing to the familiar formulae (1) only 
once the harmonic condition is reintroduced. It is worth 
noting, however, that although the generalized forms of 
the tensors will be useful for technical reasons in later 
sections, they will not give us any further physical infor- 
mation than their restriction to harmonic gauge (1). This 
is because in order to interpret t^ v and s a (JV physically, 
we must first extinguish their gauge freedom; the only 
way to do this that produces sensible local properties 
(positive energy-density, causal energy-flow, and spatial 
spin) is by insisting on transverse-traceless gauge, which 
obviously ensures that the harmonic condition is satis- 
fied. 

In addition to relaxing the harmonic condition, it will 
also be convenient to ignore matter (T^ = 0) and work 
with gravity in vacuo for the entirety of this section. 
Even though the framework of our previous papers was 
developed around the exchange of energy-momentum and 
angular momentum between matter and gravity, here we 
will be able to construct t m „ and s a ^ v from the dynamics 
of the gravitational field alone. 

A. The Fierz-Pauli Lagrangian 

We begin in a flat background spacetime (M., g a b) with 
the Fierz-Pauli Lagrangian for a massless spin-2 field [7] : 

£ FP = -L (d^Kpd^ - 2d^ a d v h\) . (3) 

From a non-gravitational standpoint, this Lagrangian 
can be derived by demanding invariance under the mass- 
less spin-2 gauge transformation [8]: 

Sh^ = => 5£fp = surface terms. (4) 



The primary aim of this section is to reproduce the 
formulae obtained in [1, 2], 

= jdf,h a/3 d u h al3 , (la) 
KS a flu = 2h 0[u dH„/\ (lb) 



2 It will often be convenient to express the formula for t^i, in terms 
of its trace-reverse r M „, as we have in equation (la). Although 
the algebraically simpler object r M ,y contains all the structure 
needed to reconstruct the energy- momentum tensor t^, it does 
not itself carry a physical interpretation. 
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For our purposes, however, it suffices to observe that Cpp 
correctly reproduces the linearized vacuum Einstein field 
equations: 



= 



SCfP l G^h a p, 



(5) 



where 

G^Kp = G« = d a d ( ^h u) a - d 2 h^/2 - d^h/2 

+ VfiV (d 2 h-d a d fS h af) )/2 (6) 

is the linear part of the Einstein tensor when the physical 
metric g ao is perturbed according to 



4>*9ab = ^ab + h a b- 



(7) 



As before, <j> '■ M. — > M maps the physical spacetime 
(■M-,g a b) to the background. 

We wish to obtain an energy-momentum tensor and 
a spin tensor from Cpp by "gauging" the symmetries of 
the background spacetime, and taking functional deriva- 
tives with respect to the background gauge fields. In the 
standard formulation of general relativity, the rotational 
and translational symmetries of Minkowski spacetime are 
gauged into a single local symmetry: diffeomorphism 
gauge invariance. As a result, the energy-momentum and 
angular momentum of matter are all contained within a 
single Belinfante energy-momentum tensor T^ el (written 
simply as T Q f, in [1, 2]) which one derives according to 
Hilbert's definition:!^ 1 = (l/^g~)(5£ maUcl /5g ab ). In 
order to obtain two separate tensors, t^ v and s°^, fol- 
lowing a Hilbert-like approach, one must gauge the rota- 
tional and translational symmetries separately^ and take 
derivatives with respect to the two gauge fields (A13). 
These are the techniques of Einstein-Cartan gravity (as 
formulated by Kibble and Sciama [9-15]) which we sum- 
marize in appendix A. 

To apply these techniques to the problem at hand, 
we will first need to covariantize Cpp. Invoking a back- 
ground tetrad and spin connection (i/", we write the 
Fierz-Pauli Lagrangian in terms of quantities which are 
covariant under local translations and rotations: 

£ FP = ± (/)„/,.,,./';"/,"• - 2D p h^D p h\) , (8) 

where D a is a covariant derivative with connection u>J Ly 1 
the volume element e = 1/ det(e^), and Greek indices 
with overbars enumerate the components of tensors in 
the non-holonomic basis {e?}. 

In order to perform arbitrary infinitesimal variations 
in and we will inevitably explore backgrounds 

with non-zero curvature R a bp,v and torsion T a p.p. For this 
reason, we must also decide how our Lagrangian should 
change when the background is no longer flat and torsion- 
free. 3 The obvious response to this uncertainty is to fol- 
low the "minimal coupling" maxim, and insist that the 



Lagrangian remain as it is in equation (8) even when the 
background is curved and contorted. Despite the sim- 
plicity of this approach, the Lagrangian Cpp is actually 
a highly unnatural choice, as it deprives the field theory 
of its spin-2 gauge-invariance when the background is no 
longer flat. To see this, consider a curved vacuum back- 
ground (Rap = 0, T a p,a = 0, Rabfiv ^ 0) and perform a 
covariantizcd spin-2 gauge transformation: 



5hp,iy = D {p £ e 



SC' F p — t) 11 ^ R a n Rh a ^ 
k ^ 

+ surface terms. 



(9) 



Thus, C' F p loses its spin-2 gauge invariance when one 
tries to extend the theory "minimally" beyond the flat 
background. 

The gauge invariance of the field theory can be pre- 
served, for vacuum backgrounds at least, if we allow h^ v 
to couple directly to the curvature of the background. 4 
The Lagrangian (10) we will use to generate t m „ and s a „„ 
will do exactly this, although we should explain that it 
is not unique in this regard. If we had wanted to present 
the subsequent calculation as a genuine ab initio deriva- 
tion of t^ u and s a ^ v , then we would need to justify our 
specialisation to (10) over the other possibilities. How- 
ever, we have already derived t^u and s a ^ from more 
concrete considerations [1, 2], and our aim here is only 
to show that a Lagrangian exists from which the tensors 
can be obtained. We will explore this curvature-coupling 
freedom in section II C, and by the end of section IV B we 
will be in a position to look back at C and better under- 
stand the significance of our "choice" . For now, we shall 
simply write down our Lagrangian as an ansatz, justified 
by its being a covariantization of £fp which preserves 
the field theory's gauge-invariance beyond the fiat back- 



the variations in and U]^" be constrained to backgrounds with 

Rabfiv 



T a ^ u = 0. However, one must enforce this restriction 



by including Lagrange multiplier terms \ R a bp,v + ^a^" T a ^ v 
in the Lagrangian, which inevitably produce unconstrained su- 
perpotentials proportional to dpX^v^ in the energy- momentum 



tensor and dp\ 



A 



in the spin tensor. Thus, one can- 



3 One might hope to avoid this decision by instead demanding that 



not use these "constrained variations" alone to construct r Ml/ and 
s0l uv fr° m a Lagrangian, as one must make further specifications 
of the form A = . . . before the energy-momentum tensor and spin 
tensor arc well-defined. In contrast, "unconstrained variations" 
can be used to construct and s a M „, as we will soon see. The 
key to this approach is that we are forced to describe how C' F p 
behaves when the background is curved and contorted. Because 
there is a precise correspondence between background-coupling 
and supcrpotcntials (see section II C) once this behaviour of the 
Lagrangian is fixed, it is able to define the energy-momentum 
tensor and spin tensor unambiguously. 
4 No set of curvature-coupling terms (nor torsion-coupling terms) 
can extend the theory's gauge invariance to include non-vacuum 
backgrounds. This comes as no surprise, considering that we 
are studying a Lagrangian Cpp that does not include matter. 
Evidently, the linearized vacuum field equations (5) can only be 
expected to be consistent when they describe perturbations from 
a vacuum background. 
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ground, and proceed to calculate its energy-momentum 
tensor and spin tensor. 

B. Energy-Momentum Tensor and Spin Tensor 

Let us consider the following the Lagrangian for the 
linearized gravitational field: 

+ •-'/'"•"/,'„„,,,/'"') . (10) 

This clearly reduces to the Fierz-Pauli Lagrangian (3) 
when the background is flat and torsion-free, and further- 
more, successfully extends the spin-2 gauge-invariance of 
the theory to curved (vacuum) backgrounds: 

Sh^ = D^v) =>■ SC = surface terms. (11) 

Treating the fields {hP" , e?, as independent vari- 
ables, we shall evaluate the energy-momentum tensor 
and spin tensor of C according to their definitions from 
Einstein-Cartan gravity: 

^Gi^)-*' sV ^Gto7^)^' (i2) 

where the subscripts e = 6 and Cj = signify that, once 
the functional derivatives have been taken, we evaluate 
the tensors on a flat torsion-free background, and the 
tetrad and spin connection become trivial (A15) to re- 
flect this. Substituting (10) into (12) we arrive at the 
following formulae for the energy-momentum tensor and 
spin tensor of the linearized gravitational field: 

nf^ = \dfj,h a ^dji ap - ^dfiKudph 013 , (13a) 

= 2h 0[u d^\f + sghrfajPp. (i3b) 

This is precisely the result we needed: C has generated an 
energy-momentum tensor t^ v and spin tensor s a „„ which 
reduce to the familiar formulae (1) when the harmonic 
condition (2) is reintroduced. 

We have achieved the main aim of this section, demon- 
strating that our energy-momentum tensor and spin 
tensor can be identified as translational and rotational 
Nocther currents of a Lagrangian for linear gravity. In 
addition, the equations (13) reveal how our tensors (1) 
generalize beyond harmonic gauge. Before we study these 
generalized tensors in detail, we shall first examine the 
freedom that was available in our choice of covarianti- 
zation of £fp, and demonstrate that the formulae (13) 
constitute a suitably unique extension of (1). 

C. Background Coupling and Superpotentials 

We begin by considering the most general Lagrangian, 
quadratic in h^ v and second-order in derivatives, which 



differs from the minimally coupled Lagrangian (8) only 
by terms which couple to background curvature; ig- 
noring surface terms, this is 

C^CJ^ + ^Rj 9 ^-^ (14) 

where £«V = = \* = is a local 

quadratic Lorentz-covariant function of h^vi the general 
form of which can be parametrized by five dimensionlcss 
constants {A„}: 

+ S^ ] (A 4 h 2 + A 5 h lS h~' s ). (15) 

If we recall the behaviour of £' FP under a spin-2 gauge 
transformation (9), it is immediately clear that the field 
theory will retain its gauge invariance for curved vacuum 
backgrounds (11) if and only if A\ = — 1. 

Inserting Cfj into (12), we find that the energy- 
momentum tensor of this Lagrangian is identical to the 
tensor (13a) derived from £, but that the spin tensor is 
given by 

+ d /3 £ a V- (16) 
In harmonic gauge this becomes 

Ks a flv = h^d a h/ + d^ lw , (17) 

revealing that {A n } are the very same constants that pa- 
rameterized the superpotential freedom of s a „„ in section 
III of our previous paper [2]. There, the value A\ = — 1 
was fixed by demanding that s a av = for all transverse- 
traceless h^ u . Now we see that this special value of A\ 
has yet another significance: it ensures that the spin-2 
gauge-invariance of the linearized theory extends beyond 
the flat background. 

Equation (17) also demonstrates that the parameters 
{A n } must take the values 

At = -1, A 2 = 1, A 4 = -1/4, A 3 = A 5 = 0, (18) 

(as they do in [2]) if the spin tensor (16) is to reduce to its 
original form (lb) in harmonic gauge; thus the freedom 
to add curvature terms (14) cannot, by itself, produce 
any other generalization of r Mt/ and s a ^ v than (13). 

Now that we understand the role played by curvature 
terms in the Lagrangian, we must also explore the possi- 
bility of coupling h^u to background torsion. If the La- 
grangian is to remain quadratic in h^ v and second-order 
in derivatives, the only contribution we need to consider 
is 

A£^--T^£/ P , (19) 



5 



where YiJ 11 ' = — ^ a vf " is composed of terms of the form 
hDh. The torsion terms generate the superpotential free- 
dom of the energy-momentum tensor: 

K^=(^ 5 J£) = W, (20) 
V2e def, J e=8 

(2=0 

which is also accompanied by a change in the spin tensor, 

" as -=ttl^)«= 2EM "' (21) 

LJ — O 

Because the energy-momentum superpotentials are of the 
form d(hdh), containing tensors of the form hd 2 h, their 
addition has the potential to spoil the homogeneous dif- 
ferential structure of (13a): r^u ~ dhdh. In fact, there 
is no superpotential e?„X '" / , entirely composed of terms 
dhdh, which vanishes in harmonic gauge. [To prove this, 
construct the most general 2-index Lorentz-covariant ten- 
sor, composed entirely of terms of the form dhdh, which 
is at least linear in d^h^ , and suppose that it is also a 
superpotential: kAt^ = dah a ^{G\dph v ji + Cjd^h^p + 
C^dph^v + T),u,(Cidph + C^d^h 1 p)) + C§d a h a ^dvh + 

Cidah a yd^h + C$d a h a ^dph? v , where {C„} are arbitrary 
dimcnsionlcss constants. Equation (20) informs us that 
d^Ar^" = for all h^; the only values of {C„} consis- 
tent with this are C n = 0.] Thus, the freedom to add tor- 
sion terms to the Lagrangian is nullified by our insistence 
that the generalized t^ v be free of second derivatives, and 
reduce to our original formula when the harmonic condi- 
tion is enforced. 

We therefore conclude that our Lagrangian (10) is the 
unique covariantization of £fp, quadratic in h^ v and 
second-order in derivatives, which according to the defi- 
nitions (12) generates an energy- momentum tensor that 
is free from second derivatives, and an energy-momentum 
tensor and spin tensor which agree with our original for- 
mulae (1) in harmonic gauge. Consequently, the result- 
ing energy-momentum tensor and spin tensor (13) are 
the unique extensions of t^ u and s a „„ beyond harmonic 
gauge, which can be derived from a covariantized Ficrz- 
Pauli Lagrangian according (12), and which do not in- 
troduce terms of the form hd 2 h into t^ v . 

Having demonstrated that (13) are indeed the unique 
extension of (1) beyond the harmonic gauge, it will be 
useful to construct their Belinfante tensor. 



D. Belinfante Tensor 

As an alternative to our description of gravitational 
energy-momentum and spin in terms of two separate ten- 
sors, t m „ and s a we may construct a Belinfante energy- 
momentum tensor [16], 

t/j,v = I'm, + d a (s fll/ + — s ^„)/2, (22) 



which combines the two. This tensor is symmetric by 
virtue of the field equations, 

*[H = T \M - 9 a s a flu /2 

= -ha^G^hfr = 0, (23) 
and is also conserved: 

0*V = = -^(W aP )G a g' r %s = 0. (24) 

Furthermore, provided surface terms are negligible, the 
Belinfante tensor defines precisely the same global mea- 
sure of energy, momentum, and angular momentum as 
t mv and s"^: 

J t ll °d 3 x = J r/d 3 *, (25a) 
J 2x [tI U ] °d 3 x = J (2x [fl r v] ° + sV)d 3 x. (25b) 

The advantage of this Belinfante description is obvi- 
ous: it combines energy-momentum and spin into a single 
symmetric tensor. 5 This apparent simplicity comes at a 
high price, however, because although the global picture 
remains intact (25) the Belinfante tensor is unable to re- 
produce the physically sensible local description that 
and s a „„ provide. 

In general, the intermixture of spin and energy- 
momentum in (22) prevents us from localizing the two 
quantities separately, and we are left with angular mo- 
mentum currents x^t^ a which "contain" spin but do 
not assign it a local current, and energy-momentum cur- 
rents i M Q which display negative energy-densities and 
non-causal energy-flux. Furthermore, because the grav- 
itational Belinfante tensor has no special geometric or 
algebraic properties in either harmonic or transverse- 
tracclcss gauge, it becomes impossible to justify a nat- 
ural gauge-fixing program. The tensor can then be 
evaluated over the entire gauge space of h^ v , and will de- 
pend on the arbitrary mapping : M — > M. as much as 
it depends on the physical properties of the gravitational 
field. 

For these reasons, we cannot advocate interpreting 
as the "true" energy- momentum of the gravitational field. 
The tensors t^ v and s a are the local measures of gravi- 
tational energy-momentum and spin, describing positive 
energy-density, causal energy-flux, and spatial spin; i M „ 
is a derived quantity which packages spin and energy- 
momentum into a single object, losing some local infor- 
mation in the process. The main application of the grav- 
itational Belinfante tensor will arise in the next section, 



5 That said, as far as our framework is concerned, the symmetry of 
the Belinfante tensor is not particularly impressive: our energy- 
momentum tensor t m „ is already symmetric, by virtue of the 
harmonic condition (2) rather than the field equations. 
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where we will also gain some insight into its physical 
interpretation. In brief, we will see that f jn , appears as 
the quadratic contribution to the Einstein field equations, 
generating perturbations in the metric alongside the (Be- 
linfante) energy-momentum of matter. In other words, it 
is the particular combination of energy-momentum and 
spin, + d a (s^ a + s^ Q - s a ^)/2, that curves phys- 
ical spacetime in a quadratic approximation to general 
relativity. It would be implausible to expect alone 
to fulfill this role, as there is no other field equation in 
which s^,, could act as the source; only by considering 
perturbations in the Einstein-Cartan equations, as we do 
in section IIIB, will we find a setting in which and 
s a '■ arise as the self-interaction source-terms in two sep- 
arate field equations. 

Proceeding with the calculation, we substitute the gen- 
eralized tensors (13) into the definition (22) in order to 
obtain the Bclinfantc tensor associated with our frame- 
work. The resulting formula can be expressed most com- 
pactly in terms of the trace- reverse of t^: 

+ \d a h p ^ {d?Kf - d v) h a P) 

+ h p{ ^d a K) a + h a[lJ G^ aPl h Pl . (26) 

Although the last term can be removed by applying the 
field equations, we will retain it for the sake of generality. 

This concludes our analysis of the Lagrangian formu- 
lation of t^v and s a ^ v . Armed with the results of this 
section, we are now in a position to "weigh" the gravita- 
tional field, and investigate the role our tensors play in 
the non-linear field equations. 



III. SELF-INTERACTION IN THE 
GRAVITATIONAL FIELD EQUATIONS 

Here we examine how t^ v and s a „ v (in their general- 
ized form (13)) occur as the quadratic terms in a pertur- 
bative expansion of the Einstein field equations and also 
the Einstein-Cartan field equations, generating gravity in 
exactly the same fashion as material energy-momentum 
and spin. 

As we move from a linear theory of gravity to a 
quadratic one, it will become important to fix the def- 
inition of h^y more precisely. Until this point, has 
been used to signify a perturbation in the physical met- 
ric: 

<fi*9ab = 9ab + Kb- (27) 

However, because our framework is based exclusively on 
gravity in the linear approximation, we could have de- 



fined h^y such that 

<\>*9ab = gab + Kb + 0(h 2 ), (28) 

and arrived at the very same results. For instance, sup- 
pose we had decided to work with the field h' v that de- 
fines a (negative) perturbation in the inverse metric: 

0*5 ab = 9 ab - h' ab . (29) 

The equation for the metric would then have been 

<$>*9ab = 9ab + h' ab + h' ac h\ + 0(h 3 ), (30) 

instead of (27), but because h' = h^ u + 0(h 2 ), the lin- 
earized theory of ft.' would be the same as h^, and our 
framework would assign the same tensors (13) to describe 
its energy-momentum and spin. Only once we came to 
study the quadratic approximation of the field equations, 
as we do now, would any mathematical difference be- 
tween the fields and h' have been observed. 

For the sake of concreteness, we will start with the 
standard definition of the gravitational field (27) and con- 
sider this to be true to all orders of approximation. As we 
will soon see, however, it is precisely the freedom to make 
field redefinitions of the form (28) that will allow us to 
cast T^ v and s a ^ v as the sources of the gravitational field; 
by the end of the next section, we will have uncovered 
a new definition of valid to quadratic order, that is 
specially selected by our local description of gravitational 
energetics. We will examine the wider significance of this 
definition, and the effects of field redefinition in general, 
in section IV. 



A. The Einstein Equations 

Consider the vacuum Einstein field equations (in the 
physical spacetime) expressed in terms of the "mixed" 
Einstein tensor density: 

V=9G a b = 0. (31) 

Mapping this equation to the background, we apply (27) 
to every instance of (j)*9ab, and expand the result in pow- 
ers of h^v : 

G^ aP Kp + G$ + 0(h 3 ) = 0, (32) 

where 

G^[r(V^O] (2) (33) 

is the quadratic part of the mixed Einstein tensor density. 

Let us now redefine the gravitational field h^ v by mak- 
ing the replacement 

V -> V + K a h a u /2; (34) 
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this gives rise to a corresponding change in the definition 
of the metric, 

(f>*9ab = 9ab + Kb + h ac h c b /2, (35) 

and causes the vacuum equations (32) to become 

G^Kp + G^{h ai h\)/2 + G$ = 0, (36) 

when working to second order. 6 Moving all quadratic 
terms to the right-hand side, and making use of the fol- 
lowing identity, 

-G<$ - G^(h aj h\)/2 = Kt^, (37) 

which is derived in appendix B, we find that the quadratic 
vacuum field equations (36) are equivalent to 

where i M „ is the Belinfante tensor of the gravitational 
field (26) constructed from and s a M „. 

Equation (38) is exactly what we had hoped to find: 
gravitational energy-momentum generates gravity in ex- 
actly the same fashion as material energy-momentum. To 
make this comparison transparent, we remind ourselves 
of the non-vacuum field equations at linear order: 

G^hap = kT^\ (39) 

where T^ 1 (written simply as in [1, 2]) is the Be- 
linfante energy-momentum tensor of matter, mapped to 
the background. Because T^ 1 is Belinfante, any intrin- 
sic spin carried by matter must be packaged inside this 
tensor according to the same formula (22) that defines 
the Belinfante tensor of the gravitational field, making 
the analogy with i M „ extremely close. Furthermore, if 
one assumes that T^ 1 is of the same order of magnitude 
as i M „ ~ 0(h 2 ), then at quadratic order the non- vacuum 
version of (38) is in fact 

d lu , a f , h aP = K{t lu ,+T^), (40) 

wherein the source of the gravitational field is the sum of 
the material and gravitational Belinfante tensors. 

It goes without saying that equation (38) can also be 
written as 

G^ a0 h a p = k (v, + d a (s^ a + s vlx a - sV)/2) , (41) 

making the function of t^ v and s a „„ absolutely clear: 
these tensors do not simply constitute a passive "kine- 
matical" description of gravitational energy-momentum 
and spin, they actively determine the field's dynamics. 



To clarify: the tensor Gjj,J is still the quadratic part of 
4>*(.V~9G a b ) when the metric is expanded according to (27); 
this tensor has exactly the same formula in terms of the new h^v 
as it did the old because the replacement (34) only alters G^i 
by quantities 0(/i 3 ), which we neglect. 



Despite the satisfying simplicity of this result, equa- 
tion (41) is clearly not the best point at which to end our 
investigation. Having extolled the virtues of a formal- 
ism which keeps spin separate from energy-momentum, 
our real goal must be to find a formulation of gravity in 
which Tjj, v and s a „ v appear as source-terms in separate 
gravitational field equations. It should come as no sur- 
prise that Einstein-Cartan theory will provide precisely 
the environment in which to achieve this objective. 

B. The Einstein-Cartan Equations 

We will now disentangle the spin and energy- 
momentum in equation (41), formulating a quadratic ap- 
proximation to Einstein-Cartan gravity in which t^ v and 
s ° ' uv a PP e & r as separate source-terms in the field equa- 
tions. In close analogy with the previous section, we pro- 
ceed by expanding the field equations (A12) to second 
order in /„„ and w t /" / , where 

<t>"e% =51- /y2, (42a) 
<jr U J»=w t r, (42b) 

are initially considered to be true to all orders; we then 
perform a non-linear field redefinition, 

/„„ + 0(/ 2 ), (43a) 

V" V" + °(/ 2 )' ( 43b ) 

to generate the field equations we desire. In order to 
identify and s a „ v in these equations, it will also be 
necessary to express tensors of the form wdf + w 2 in 
terms of /i M „. To this end, we will evaluate these tensors 
on torsion- free perturbations (A22), 

= (d [v r ] a + d^f a ^ + d a f^)/2 + O(f), (44) 
in the "symmetric" rotation gauge (A24): 

f\M = 0{f)- (45) 
These relations allow us to identify 

f M u = V + 0(f 2 ), (46a) 
= d^h^ a + 0{f), (46b) 

and thus convert the quadratic parts of the field equa- 
tions into the corresponding tensors of perturbative gen- 
eral relativity: wdf + w 2 = dhdh + 0(/ 3 ). 7 



7 The validity of the formulae (13) can only be guaranteed in the 
generally-relativistic regime, i.e. torsion-free gravity described by 
a symmetric tensor field h^. In the linearized Einstein-Cartan 
theory, this corresponds to the restriction (44) for •w c ! il ' and (45) 
for ftiv If we already knew how to generalize and s a M „ be- 
yond the generally-relativistic regime (where the gravitational 
field is described by unconstrained / M „ and «J a M ") then it would 
be possible to recognise these tensors in the field equations with- 
out making such restrictions. The conversion dh — > df + w is 
not unique, however, so it is not immediately clear how this gen- 
eralization should be achieved. 
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To begin, let us focus our attention on the first 
Einstein-Cartan field equation (Al2a). Following the ap- 
proach of section III A, we express the vacuum field equa- 
tion in terms of the mixed Einstein tensor density, 

eG a b = 0, (47) 

where G a b = (Rj L —e l R/2)e b ll is a function of the physical 
tetrad and physical spin connection wj 1 " . Mapping 
this equation to the background, we expand to quadratic 
order in f^ v and wj 11 ' , and simplify the resulting equa- 
tion by taking the trace-reverse: 

+ f a fS d [li w a]v P-2w M a v w m P a . (48) 
We now redefine w a ^ v according to the replacement 

v *>cr - fw a ^/2 + f^w a ^ 

+ f^d a f^/4; (49) 
the quadratic field equation (48) then becomes 

2d [tl w a]u a = d [fl fw a]u a - <9[ A1 |/„ (3 w N ' 9Q - d [t ,f afi w a]vf} 

-2w M % Wl / a -d M f v d la] r /4 

+ d^f a d a] .U /4. (50) 

Applying (46) to the terms on the right-hand side and 
taking the trace-reverse of the equation, we obtain 

2d[i t w a ] l , a - rj^vdaWp 010 = kt^. (51) 

This is the field equation we had hoped to construct, mir- 
roring the structure of the linearized non-vacuum field 
equation (A21a) with gravitational energy-momentum 
t m „ taking the place of the material energy-momentum 
tensor T^ v . 

We now turn to the second vacuum Einstein-Cartan 
field equation (Al2b). Writing this as 

T% = (52) 

in the physical spacetime, we once again use the stan- 
dard definitions (42) and expand to second order in the 
background: 

d [fl f a v] + 2w [ ,% ] = f^r^/2 + w p a {v ft A 

+ WV (53) 

Consistency with the first field equation (51) requires us 
to redefine wj 11 ' as before (49) but places no constraint 
on the definition of ; we are therefore free to make the 
replacement 

Uu / M „ - WV4- (54) 

Appling these redefinitions to the second field equation 
(53) and converting the quadratic terms using (46), one 
finds that 

d^r*] + = 2h p[u d^h^ 



from which the desired equation follows: 

d [tl r v] + 2 W[ « V] 

+ ^\(9W)f - dpfw] ~ 2w/ M ) = ( 56 ) 

We have found a suitable counterpart to equation (51), 
in which the gravitational spin tensor s a ^ takes on the 
role played by material spin S a „ u in the linearized field 
equation (A21b). 

Combining these results, we conclude that under the 
perturbative expansions 

</>*e?=^-/y2 + rjV8, (57a) 

+ fp^dafV/A, (57b) 

the vacuum Einstein-Cartan field equations are approxi- 
mated, to quadratic order, by 

2d[^w a ] u a - ■q tiu d a w f} al3 = kt m „, (58a) 

d {li r v] + 2 W[ « V] 

wherein t^ v and s a generate the gravitational fields 
f^v and wj 11 ' in an identical fashion to the energy- 
momentum and spin of matter (A21). We have found the 
analogue of (38) in the Einstein-Cartan theory of gravity, 
in which gravitational energy-momentum and spin act as 
the source terms of separate field equations. 

IV. FIELD REDEFINITION 

We have succeeded in demonstrating that t^ v and s a ^ v 
do indeed express the dynamical "weight" of the gravi- 
tational field, and have uncovered the field definitions, 
(35) and (57), which make this relationship manifest at 
the level of the field equations. We now turn our atten- 
tion to the field definitions themselves, investigating the 
importance of (35) and (57) in a broader context, and 
exploring the effects of field redefinitions in general. In 
particular, the analysis of section IV C will allow us to 
confirm the statement that justified the traceless condi- 
tion, which we imposed when deriving the formula for 
s a ^ u in our previous paper [2]. 

A. The "Central" Expansion 

In many respects, the most striking aspect of the new 
definition of h^, as displayed in (35), is how closely it 
relates to a linear perturbation in the metric (27) and a 
linear perturbation in the inverse metric (30); this is in 
comparison with the full range of local Lorentz-covariant 
field definitions consistent with (27) at linear order: 

<t>*9ab = 9ab + h a b + Bih ac h c b + B 2 h ab h 

+ 9ab(B 3 h 2 + B 4 h cd h dc ) + Q{h 3 ), (59) 
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where {B n } arc arbitrary constants. One can argue, in 
fact, that the definition (35) lies at a natural "center" of 
the four-dimensional space parameterized by {£?„}. This 
argument begins by observing that a priori there is no 
special variable which represents the "true" dynamical 
field of general relativity: one can equally well define the 
gravitational field as a linear perturbation in a metric 
density (—g) x g a b, or an inverse metric density (~g)' K g ab , 
for any value of A. Of all these choices, perturbations 
in the metric and its inverse (i.e. A = 0) arc distin- 
guished by the fact that they possess linear-order gauge- 
transformation of the form without a part propor- 
tional to T]^ v da^, a , and can therefore be identified with 
the Fierz-Pauli massless spin-2 field. However, once we 
have restricted our interest to these particular definitions 
((27) or (30)) the decision to focus on one, and discard 
the other, is completely arbitrary. Instead of making a 
forced choice between two essentially equivalent options, 
one might instead consider the definition that lies exactly 
half-way between them, where the values of {B n } are the 
mean of those in (27) and (30). It is easy to see that this 
"center point" is precisely the field (35) that casts as 
the self-interaction term of the quadratic field equations 
(38)! This is an extraordinary coincidence, as t^ v and 
s a were selected for their capacity to display positive 
energy-density, causal energy-flux, and spatial spin; none 
of these criteria would be expected to determine a defi- 
nition of h^y that is geometrically distinguished in this 
way. 

B. Expansion of the Einstein-Hilbert Lagrangian 

The new metric expansion (35) also offers some per- 
spective on the particular form of the Lagrangian (10) 
that generated the main results of section II, including 
the generalized formulae (13) for and s a „„. 

In an earlier paper [17] we studied the expansion of the 
Einstein-Hilbert Lagrangian, 8 

£eh = -V=gR/n, (60) 

under a linear perturbation of the inverse metric: h' 
as defined in (29). 9 Taking care to retain all terms pro- 
portional to background curvature, but ignoring surface 
terms, the Einstein-Hilbert Lagrangian was found to ex- 
pand as follows, 

£ E H[.9 afc - h lab ] = -^gR/K + £' 1 [ti ab ] + C' 2 [h' ab ] 

+ 0(h' 3 ), (61) 



As it was Hilbcrt alone who formulated general relativity in terms 
of a least-action principle [18], some authors refer to (60) as the 
Hilbert Lagrangian. 

The field written as h^y in [17] is in fact —h'^ in our present 
notation; the Lagrangians of that paper also take the opposite 
sign to those here. 



where the linear and quadratic parts of the Lagrangian, 
£[ and C 2 , are given by equations (65) and (66) of [17]. 

To generate the expansion of £eh under our newly de- 
fined perturbation h^ v , we need only compare its metric 
expansion (35) to that of h' (30), 

h' ah = h ah - h ac h c b /2 + 0{h 3 ), (62) 

and substitute this relation into (61). After commuting 
derivatives, and discarding surface terms, the quadratic 
part of this expansion turns out to be 

C 2 [h ab ] = C' 2 [h ab ] + C' 1 [-h ac h c b /2] 

= ^ {V c h' ab V c h' ab - 2V a h' ab V c h' b 

+ 2h' ab R cabd h' cd ) , (63) 

which is precisely the form of the Lagrangian C that we 
used to reproduce the formulae for and s a „„! Thus, 
the curvature term in (10), which we introduced in sec- 
tion II as an ansatz, can be understood as a consequence 
of the special definition of hn V associated with our frame- 
work: these are simply the terms proportional to the 
background curvature that appear when the Einstein- 
Hilbert action is expanded to quadratic order. 

C. Field Redefinitions and Superpotentials 

To understand the new metric expansion (35) in a 
wider context, we should also explain the relationship 
between field redefinitions and the superpotentials we en- 
countered in section II C. 

Recall that the addition of superpotentials to and 
s a „ v corresponds to the addition of curvature terms (14) 
and torsion terms (19) to the Lagrangian. Notice, how- 
ever, that is unaffected by torsion terms: the energy- 
momentum superpotentials (20) cancel those of the spin 
tensor (21) when they enter the formula (22). As a result, 
the superpotentials of the Bclinfante tensor arc charac- 
terized by curvature terms alone, which were determined 
by the five parameters {A n } of equation (15). Although 
these parameters are fixed according to (18), for the sake 
of argument let us relax these equations and alter each 
A n by an amount AA n ; the gravitational Belinfante ten- 
sor then gains the supcrpotcntial term 

«Ay = Kd a {& s ™ + A S y - A*v)/ 2 

= d Q ^A£ Q V, (64) 
where, according to (15), 

A£«V = AA^fXj + AA a W#«Sg + AA 3 hj/^ 

+ 8fX\ (&A 4 h 2 + AA 5 h lS W s ) . (65) 

It is also possible to generate superpotentials in the 
quadratic approximation to Einstein's field equations: an 
arbitrary field redefinition 

hpv hfw + Ah,i V , (66) 
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adds the divergence-free tensor 

Q/3 A/i a/ 3, 



(67) 



to the right-hand side of the field equations (38) and thus 
defines a new Bclinfante tensor, 



^UU — K tfJ,u G„ v Ah a 



(68) 



that acts as the source of the new gravitational field. 
Therefore, as long as we can find a field redefinition Ah^ v 
to solve 



GS^AhctB = d a d p AY, av 



(69) 



we can produce the same superpotential in the field equa- 
tions as the ones we have generated by altering {A„} in 
the Lagrangian. 

First we shall try to solve equation (69) using local 
field redefinitions. Noting that Ah^ v will also need to 
be Lorentz-covariant and quadratic in to solve this 
equation, the most general field redefinition we need to 
consider is 



Ah,,* = AB 1 h iia h a v + AB 2 h^h 

+ r 1 ^(AB 3 h 2 + AB i h al3 h a ?) 



(70) 



where the {AB n } correspond to changes in the param- 
eters {B n } of equation (59). Comparing the number of 
free parameters here with those of (65) it is immediately 
clear that these local redefinitions will not span the entire 
space of superpotentials. Indeed, if we use the identity 
G™PX a p = -2d a d (6tx*h valid for any symmetric 
tensor X^ v , we can rewrite (69) as 



= d a d f3 (A£ tt V - 25^A^J 



(AA 2 - 2AB 2 )hhp"J ] 



(71) 



(72) 



inserting (70) and (65), we arrive at 
= d a d f) lAAxh^h 

+ (AA 3 - 2AB 1 )hj/ffl 

+ S^ ] (AA i + 2AB 3 + AB 2 )h 2 

+ S^(AA 5 + 2AB 4 + AB^sh-' 6 



which makes the mismatch of parameters unequivocal. 
Clearly, this equation can only hold for all if 



AAi = 0, 



(73) 



and if we assume that this is the case, the local field 
redefinition Ah^ u is determined uniquely: 

ABi = AA 3 /2, AB 2 = AA 2 /2, (74) 

AB 3 = -AAi/2 - AA2/4, ABi = -AA 5 /2 - AA3/4. 

As the AA n were defined relative to the values (18) that 
correspond to our Bclinfante tensor t^ v , the condition 
(73) implies that 



Ax 



-1, 



(75) 



which also arose in section II A as the requirement that 
ensured the curvature terms would extend the spin-2 
gauge invariance of the Lagrangian beyond the flat back- 
ground. 

If A\ 7^ —1 then it is still possible to solve equation 
(69) by inverting the differential operator G va ^: 



ah; = -d Q ^A£-V, 



(76) 



where the precise form of the propagator l/<9 2 will de- 
pend on boundary conditions. As we have seen, however, 
this Ah^y cannot be a local function of h^ v ; hence it will 
no longer be possible to express the physical metric (j)*9ab 
as a local function of the gravitational field (59). 

From this vantage point, we can now appreciate an- 
other important property of our framework: the Belin- 
fante tensor t^ (constructed from r^ v and s a ) is the 
source of a gravitational field of which the metric 
is a local function (35). Of course, this was precisely 
the reason for fixing Ai = —1 (by imposing the trace- 
less condition) when we derived the formula for s a (11/ in 
our previous paper [2]. Note that the statement that we 
made in section IIIB of that paper, justifying A\ = — 1, 
has now been proven. 

This state of affairs should be contrasted with that of 
the celebrated tensor i^L f Landau and Lifshitz [3] : 



Kill = -^G ah 



1 



V c V d U*(3S a[ y ]d ) 



(77) 



The divergence on the right-hand side clearly contributes 
a term 



(78) 



at second order, corresponding to a superpotential with 
AAi = 1 => A\ = 0. Thus, there can be no local field 
redefinition that will render the Landau-Lifshitz tensor as 
the source term of the vacuum Einstein field equations, 
and furthermore, the tensor cannot be derived from a co- 
variantized Lagrangian which maintains its spin-2 gauge 
invariance beyond the flat background. 

It is rather surprising that this deficiency is not more 
widely known. In their effort to construct a gravitational 
energy-momentum tensor that was symmetric and free 
of second derivatives in all gauges, it seems that Lan- 
dau and Lifshitz were forced to include a superpotential 
(78) that would be impossible to generate in the field 
equations by a local redefinition of h^ v . The only way 
the Landau-Lifshitz tensor can be given an equal footing 
with the energy-momentum of matter, generating gravity 
alongside in the Einstein field equations, would be to 
define the gravitational field h^ v in terms of a non-local 
perturbation in the metric. 



D. Beyond Second Order 

It will be difficult to gain any further insight into the 
physical meaning of the new h^ v without first deciding 
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how the definition (35) should extend beyond quadratic 
order. As this question is intimately related to the issue 
of extending the formulae of and s a ^ v to the full non- 
linear theory, we shall postpone a thorough investigation 
of this topic for another article. For now, we only men- 
tion one particularly attractive possibility. As we have 
already explained, the new definition (35) lies on a point 
of symmetry between a linear perturbation in the metric 
(27) and a linear perturbation in the inverse metric (30); 
reflecting this, one finds that the expansion of the inverse 
metric, consistent with (35) is 

0*g ab = 9 ab - h ab + h ac h c b /2 + 0(h 3 ), (79) 

which, apart from a change in sign convention /i M „ — > 
—h^ u , is identical to the metric expansion (35) to 
quadratic order. Thus, a particularly natural extension 
of (35) would be one which preserved this symmetry ex- 
actly, so that the metric and its inverse had an identical 
expansion to all orders, except for a sign change in h^. 
This idea can be realized by viewing the tensor h a b /2 as a 
linear map and forming its exponential [e' l/ ' 2 ] a b ; a metric 
defined by 

4>*9a b = [e h / 2 ] c a9c d [e h/2 ] d b , (80) 

is then consistent with (35) to second order, and more- 
over, the associated expansion of the inverse metric, 

<t?g<* = {e- h / 2 ] a c g cd {e- h / 2 } b d , (81) 

is clearly identical to (80) apart from a change in the sign 
of h^y. 

E. New Fields for Einstein-Cartan 

While it is certainly tempting to bring our analysis 
to bear on the new field definitions (57) that arose in 
Einstein-Cartan theory, unfortunately a full discussion 
of these variables will not be possible at this time. As we 
previously explained, the right-hand sides of (58) are only 
given by the formulae (13) when we restrict ourselves to 
the torsion-free perturbations (44) and symmetric gravi- 
tational field (45) of general relativity. We are therefore 
free to alter the quadratic parts of the new field defini- 
tions (57) by terms proportional to T a ^ u and /u^i: such 
terms vanish under the aforementioned restrictions, and 
so do not interfere with our results. Until we fix the 
definitions of and s a „„ in terms of arbitrary pertur- 
bations f^x, and w a ^ y , this degeneracy will remain, and 
it will be difficult to offer a physical interpretation of the 
new field variables. We leave this generalization of r^ v 
and s a „ v , and the consequent analysis of the field vari- 
ables they define, for another paper. For now, we shall 
simply remark that the new tetrad expansion (57a) de- 
fines an inverse tetrad 

<t>*< = k + F a /2 + nrjs + o(f), (82) 



which, to second order, differs from the tetrad expansion 
only by a change in sign convention f^ u — > —f^u'- mir- 
roring the relationship between the metric expansion (35) 
and inverse metric expansion (79) for the new definition 
of h^v that arose in section III A. In fact, f^ v defines a 
metric 

$*9ab = 9ab + f(ab) + f(a\ %|6) /4 + .TaW 4 

+ 0(f), (83) 

which, if we set 

ffj,v = V + 0(f), (84) 

gives 

4>*9ab = 9ab + Kb + + 0(h 3 ), (85) 

identical to the new metric expansion (35) of section 
III A. Note that the equivalence (84) between and 
/i M „ now holds to quadratic order, whereas one would 
only expect a linear correspondence from the symmetric 
gauge (A26) condition alone. 

V. CONCLUSION 

As a local description of the energy-momentum and 
spin of the linearized gravitational field, t^ v and s a ^ u 
serve a number of purposes within the theory. In addi- 
tion to accounting for the energy-momentum and angular 
momentum exchanged locally with matter [1, 2], we can 
now confirm their status as Noether currents of transla- 
tional and rotational symmetry, and as sources of gravity 
itself. Thus, our framework displays many of the funda- 
mental properties possessed by previous treatments of lo- 
cal gravitational energetics [3-5] , all the while endowing 
linear gravity with positive energy-density, causal energy- 
flux, and spatial spin. 

In summary, the main results of this article are as fol- 
lows: 

• A Lagrangian (10) has been constructed from which 
our gravitational energy-momentum tensor t^ v and 
spin tensor s a „„ arise according to standard varia- 
tional definitions (12) once harmonic gauge (2) has 
been enforced. 

• The formulae for t m „ and s a have been gener- 
alised beyond harmonic gauge (13). This is the 
unique generalisation that does not add terms of 
the form hd 2 h to t^ v . 

• The Belinfante tensor associated with our descrip- 
tion tpv = r^+d a (s^ a + s Vfl a - s a fJil/ )/2 has been 
calculated (26). 

• A non-linear perturbation in the metric has been 
uncovered (35) which generates t^ v and s a (com- 
bined into t^v) as the second-order contribution to 
the Einstein field equations (38). 
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• Non-linear perturbations in the tetrad and spin 
connection have been found (57) which generate 
Tfj, v and s™^ separately as the second order contri- 
butions to the Einstein-Cartan field equations (58). 

• The global quantities defined by t m „ and s a ^ v have 
been examined (appendix C) and shown to agree 
with ADM energy-momentum and angular momen- 
tum under appropriate conditions. 

In addition, we have described at quadratic order: 

• The correspondence between background coupling 
in the Lagrangian and the superpotentials present 
in the gravitational energy-momentum tensor and 
spin tensor that the Lagrangian defines. (Section 
II C) 

• The correspondence between these superpotentials 
and the field redefinitions required to generate 
those particular gravitational energy-momentum 
tensors and spin tensors in the field equations. 
(Section IV C) 

Finally, as a consequence of the above analysis: 

• We have isolated the range of superpotentials 
{Ai ^ —1 in (15)) which (i) cannot be generated by 
a Lagrangian that maintains the gauge freedom of 
h a b beyond the flat background (11) and (ii) can- 
not be generated in the field equations by local field 
redefinitions (59) . The Landau-Lifshitz tensor con- 
tains a superpotential of this form; as such, non- 
local field redefinitions are required to render the 
Landau-Lifshitz tensor as a source for the gravita- 
tional field. 

Having embedded t^ v and s a „ v within various aspects 
of the linear and quadratic approximations to gravity, 
the key goal that remains is to extend these ideas to the 
full non-linear theory. This is obviously an ambitious 
task, and at the present stage it is far from clear which 
of our framework's properties can survive in the exact 
theory. However, based on the results of this article, it 
seems more than likely that Einstein-Cartan gravity will 
provide a natural starting point from which to launch 
this undertaking, with the field definitions (35), (80) and 
(57) offering clues as to the new field variables into which 
this theory should be cast. 
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Appendix A: Einstein-Cartan Theory 

The role of this appendix is to briefly introduce 
Einstein-Cartan gravity, establish notation, and serve as 
a reference for results needed in the body of the paper. 
For a more complete treatment of the subject, see [12- 
15]. 



1. Kinematics 

Einstein-Cartan theory is a slight extension of gen- 
eral relativity, in which (as formulated by Kibble [9] and 
Sciama [10, 11]) translational and rotational symmetries 
are gauged separately, rather than being subsumed into 
a single diffeomorphism gauge transformation. The grav- 
itational field is represented by four vector fields (the 
tetrad) and six covector fields loJ 11 ' = ~^ a u ^ (the spin 
connection). Following the conventions of our previous 
papers [1, 2] and Wald [19], Greek letters are used as nu- 
merical indices (running from to 3, raised and lowered 
by "q^u = diag(— 1, 1,1,1)) while Roman indices represent 
the tensor "slots" of Penrose's abstract index notation 
[19, chap. 2.4]. Note that the numerical indices now come 
in two varieties: the unadorned Greek letters are used 
to enumerate the components of tensors in a Lorentzian 
coordinate system {x^} of fiat spacctime, whereas the 
Greek letters with overbars enumerate the components 
of tensors with respect to the non-holonomic basis {e?} 
formed from the tetrad. For example, 

vn=e%v a , v* = e£v a , (Al) 

where e% is the inverse tetrad, defined by e^e^ = 5£. In 
order that this convention be consistent with general rel- 
ativity (in which abstract indices are raised and lowered 
by the metric g a b) we identify the metric as follows: 

gab = e£e A5 . (A2) 

The tetrad and the spin connection are gauge fields: 
they allow the global translational and rotational symme- 
tries of fiat space to be generalized to local symmetries. 
Local translations are brought about by diffeomorphisms 
ip : M. — > Ai, under which a scalar field ip has a particu- 
larly simple transformation law: %p(x) — > ip(tp (x)). The 
tetrad allows us to form a "translation covariant" deriva- 
tive dp = e^da] as this derivative carries no spacetime 
indices, dpi/j will also transform as a scalar field. Local 
rotations are embodied by position-dependent Lorentz 
transformations which act on the tetrad, 

e£ eSA%(aO, A%A^ = r,^. (A3) 

Neither d a nor dp transform covariantly under local ro- 
tations (d a vji d a (vpA p ji ) = A v pd a Vy + vvdalsPp) so a 
"rotation covariant" derivative D a is constructed using 
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the spin connection: 



3. Dynamics 



DaVp, = d a Vfi - W^pVa, 

D a v' l = daV ii + u a ' l s v a , etc. 



(A4) 



If we declare that the spin connection should transform 
according to 



(A5) 



under local rotations, then it is easy to show that these 
derivatives arc indeed covariant: 



DaVjl 

D a v* 



(A-^W, etc. 



(A6) 



Using the tetrad once again, we can now construct a 
derivative Dp = epD a that is covariant under both local 
translations and rotations. Thus, by replacing all par- 
tial derivatives with covariant derivatives, and all space- 
time indices with basis indices, we can "gauge" the global 
Poincare-invariance of any flat-space field theory, and in 
doing so, extend the theory to a spacetimc with curvature 
and torsion. 



2. Curvature and Torsion 

The rotation covariant derivatives define a curvature 
tensor Rj fi : 

[D a ,D b ]v, = -Rj- Vfl Vt; p , (A7) 
from which it follows that 

R ab % = 2 (0[aWfc] A * + W[a|V>|6]%) ■ (A8) 

Contracting this tensor with the tetrad yields an asym- 
metric Ricci tensor R bt , = epi? ab A1 p , and Ricci scalar 

R = e-R P- 

The major difference between Einstein-Cartan gravity 
and general relativity is that, in addition to curvature, 
spacetime may also possess torsion, which we quantify 
with the torsion tensor T a ab = 2D[ a e^ , or equivalently, 



' 11.1 



(A9) 



If this tensor vanishes everywhere, we can use D^e^ = 
to relate the spin connection to the tetrad: 



HV _ Mrr vb 

2 — t- b v a e , 



(A10) 



The dynamics of the tetrad, spin connection, and 
matter fields are determined by a Lagrangian Ceg that 
closely resembles that of the Einstein-Hilbert action: 



£ec = —eR/n + £ r 



(All) 



where £ ma tter is the (covariantized) matter Lagrangian 
and e = det(e^) = 1/ dct(e^) = sf^g is the volume ele- 
ment. Variation with respect to and ooJ Xv generates 
the gravitational field equations, 



G n 



(Al2a) 
(Al2b) 



where matter's energy-momentum tensor T a ^ , and spin 
tensor iS°np, are the conjugate currents of the transla- 
tional and rotational gauge fields, 



1 SC B 



2e Se" 



cia 1 ^Anatter 



and we have written 



Rf - e£R/2, 

a b n 



F% 2e- 1 D b (eef B 4 



(A14a) 
2ef A n ]a . (Al4b) 



Consequently, the energy-momentum of matter gener- 
ates curvature, and the intrinsic spin of matter generates 
torsion. When S a n P = everywhere, the second field 
equation (A12b) ensures that torsion will vanish also; on 
substitution of (A10) and (A2), the first field equation 
(A12a) then becomes the usual Einstein field equations. 



4. Perturbations 

When the curvature and torsion of the physical space- 
time are small, it is often convenient to represent 
the gravitational fields as perturbations from a flat 
torsion-free background (M, , &J 1 "). This spacetime 
is equipped with a constant tetrad 



el = Si = (da*) OI 



(Al5a) 



defined by a Lorentzian coordinate system {x^}, and a 
vanishing spin connection, 



0. 



(A15b) 



where V a is the familiar (torsion-free) metric-compatible 
derivative of general relativity. Thus, when there is no 
torsion, the covariant derivatives Dp are equivalent to V a 
in the sense that Dpv 9 = epe^\7 a v b . Thus, in the absence 
of torsion, the curvature tensor R ab "p is equivalent to its 
counterpart from general relativity. 



In the background spacetime it is customary to ma- 
nipulate indices using the background tetrad; because 
e% = Stf, the distinction between barred indices and un- 
barred indices can then be dropped. 

Mapping the physical spacetime (M., e^, w/") onto 
the background with a diffcomorphism (j> : M. M., we 
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define an asymmetric tensor field /" as a perturbation 
in the physical tetrad, 



4>*el = S%- /y2, 



(A16) 



and the tensor field = — as a perturbation in 

the physical spin connection, 



(A17) 



The perturbation in the tetrad (A16) will be accompa- 
nied by a perturbation in the inverse tetrad, 

4>*e» = Si + FJ2 + rj\/A + 0(f), (A18) 

which in turn perturbs the physical metric (A2): 

<f>*9ab = 9ab + f(ab) + f {af f »\b) / % + f^aUb/^ 

+ 0(f). (A19) 

Thus, in the linear approximation, we can identify the 
symmetric part of / M „ with the metric perturbation h^ v 
of general relativity: 



f( P ,)+0(f). 



(A20) 



Working to first order in / M „ and w a ^ v , the Einstein- 
Cartan field equations (A12) take the following form: 

2d[ t ,w a]u a - v^daWp" 13 = kT^, (A21a) 

+ S^(d M f-d ^ M -2w/ M ) = K S% v , (A21b) 

where, for the sake of notational brevity, we have dropped 
the <)>* from the tensors 0*1^ and <j)*S a ^. To recover 
the linearized field equations of general relativity, we need 
only set S a ^ v — 0: the solution to equation (A21b) is 
then 

w/" - (d [ T ] a + d^f a ^ + d a f™)/2, (A22) 
which can be substituted into (A21a) to retrieve 

Gfj, u af " /(a/9) = kT^. (A23) 

The equivalence between / M „ and can be strength- 
ened further by fixing the rotation gauge freedom (to 
linear order) with the condition 



/[ H = 0(/ 2 ). 



(A24) 



This can always be achieved by a local rotation (A3) of 
the form 



A^ = ^ + (/'W I /)/4 + 0(/ 2 ), 



(A25) 



which is indeed a valid Lorentz transformation, 
^"(1^^10/3 = i)iiv + 0(f 2 ), and has the desired effect: 
ftiu ffrv) + 0{f 2 ). In this "symmetric" gauge, / M „ 
and h^u are equivalent to linear order: 



h fi » = Uv + 0(f). 



(A26) 



Appendix B: An Identity 

Here we derive an identity that relates the Belinfantc 
tensor (26) to the Einstein tensor of physical spacctimc. 

First, for notational purposes, let us define a tensor & 2 ^ a b 
to represent the quadratic part of the "mixed" Einstein 
tensor density: 

9G a b )f } . 



G^ b = y 
This definition expands to give 



(Bl) 



G$ = R$ R$h a v - v,AR {2) - R^h a0 )/2, (B2) 



where 



R 



(1) 



b*R ab ] 



(1) 



R 



(2) 



*Rabf\ (B3) 



are the linear and quadratic parts of the Ricci tensor, 
when expanded according to (27): 

R$ = d a d ( »K) a ~ d 2 h^/2 - d^h/2, (B4a) 

R$ = -h aP {2d a d (li h v)p - dpdrhep - dedph,,*)/! 

+ d^d„h a p/A + a a h\(d [a h p]v ) 

- d a h a P(d (ll h v)f3 - dph^/2). (B4b) 
Now consider the tensor 



(B5) 



and take its trace-reverse: 

On* = R$ - R { ^ a u + \d a d^(h v) ^) 

- ^d 2 {h m h a v ) - ^a^heph^). (B6) 

Substituting equations (B4) into (B6), one finds 
- ]^ h a ( li {d a d f3 h v) i 3 - d v )d a h) 

+ h(2d a d^K )a - d 2 h^ - a^h) 



-d^d u h a p - -d a h\d p h ua 



-d a h afs (d {ll h u)f3 - dphpv) + -d a h p[tl d v) h 
+■ - A d a h(2d { ^K )a - d a V) - h\R$ a 



a/3 



4 



(B7) 



the last line of which can be confirmed by expanding 
out all the trace-reversed fields on the right-hand side 
of (26) and observing that h a [v R^) = h a [v G [ ^ a = 

h a [ 1/ G /1 ] al3 ' y hp 7 . Inserting (B5) into the trace-reverse of 
(B7), we conclude that the following identity 

«V = "Gg ~ G^(h aj h\)/2, (B8) 
is valid for all h mJ . 
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Appendix C: ADM Energy-Momentum 

In this article, and those that have preceded it [1, 2], 
we have focussed on the local aspects of gravitational 
energy-momentum and spin; although we will not at- 
tempt a thorough investigation here, it will be valuable 
to briefly examine the global energy and momentum that 
our framework defines, and compare these quantities to 
the well-known results of Arnowitt, Deser, and Misner 
(ADM) [20-22]. 

Recall that, as seen in (25), the Belinfante tensor t^ u 
defines the same total energy, momentum, and angular 
momentum as t^ v and s™^. 10 Thus, for the purposes of 
this appendix, we are free to use whichever set of tensors 
is convenient, and the results we derive will carry over 
to the other. With this in mind, let us define the total 
energy-momentum of gravity and matter by 

= J d 3 y^(T Bc V + (r 1 )*0 (dM/T(dy%, 

(CI) 

where = (0 _1 )*a; At are the images of the Lorentzian 
coordinates in the physical spacetime. 11 To this 

definition we now apply the Einstein field equations: 



v* = -j d3 y^9 {G a b + n{r l yt a h ) ( wwdA 

= ~J d 3 x (G^Kp + G< 2 V + <) , (C2) 

where, in the second line, we have evaluated the integral 
in terms of background quantities, expanded the met- 
ric according to (27), and neglected terms 0(h 3 ) under 
the assumption that the gravitational field is everywhere 
weak enough that the quadratic approximation to gen- 
eral relativity will suffice. We now use the identity (37) 
to write the energy-momentum as 



<P„ = ~ / d 3 xG^ (h aP + 0{h 2 )) 



(C3) 



Although terms 0(h 2 ) cannot be neglected in general 
(otherwise t^ should never have appeared in the integral 
(CI) to begin with) we note that all the terms in (C3) 
are total spatial derivatives, so will depend only on 
the behavior of the gravitational field on the boundary 



10 To be precise: the integrals in (25) may in fact differ by surface 
terms quadratic in However, as we will see, these can be 
neglected in comparison to the surface terms linear in h^u when 
the boundary of the integral is taken to spatial infinity. 

11 In equation (CI), the contraction between the energy-momentum 
tensors and the covector (dy°) a defines, in the usual way, the 
energy-momentum densities on the surface of integration y° = 
const. The vectors (d/dy M ) a have assumed the role of killing 
vectors in the absence of an exact spacetime symmetry; these 
same vectors were denoted by ej^ in [I, 2], but this symbol is 
now being used for the Einstein-Cartan tetrad. 



of the integral. Thus, as the limit is taken in which this 
boundary moves to spatial infinity, we will require the 
linear surface terms dh ~ l/r 2 in order that integral be 
finite, and as a consequence the quadratic surface terms 
hdh ~ 1 /r 3 will be negligible in comparison. 

Let us first consider the total energy of the system: 

V ^ I n'.ni),;),!,,, - didihjj) 



2 k 



d Si(djhij - dihjj), 



(C4) 



which the reader will recognise as the ADM mass [21, 22]. 
Furthermore, the total linear momentum 

= — f d 3 x(d j h ij - d k d k h 0l - dji^ + didjh 0j ) 



2 k 



d 2 Sj(hij — djhoi — Sijhkk + 2dijdkhok — dihoj) 



(C5) 



which is the familiar expression for the ADM momen- 
tum [21, 22] truncated at linear order. Thus, when the 
terms 0(h 3 ) can be neglected from the field equations, 
and the terms 0(h?) can be neglected at spatial infin- 
ity, our gravitational Belinfante tensor t^ defines ex- 
actly the same total energy and momentum as ADM. 
Moreover, these results would also arise if we had de- 
fined using t^, rather than the Belinfante tensor 
tfiu. Thus, is able to cast the global information 
present in the ADM energy-momentum in terms of a local 
description with physically sensible properties, including 
gravitational energy-density that is nowhere negative and 
gravitational energy-flux that is nowhere spacelike. 

Although the ADM energy-momentum is usually rep- 
resented in terms of the asymptotic behaviour of the 
gravitational field, as above, the reader should also be 
aware that these global quantities can be cast as spa- 
tial integrals of a gravitational Belinfante tensor t^ M 
that emerges from the canonical formalism [20-22]. Al- 
though ADM did not propose that this tensor should be 
interpreted as a physically meaningful local measure of 
gravitational energy-momentum, it is nonetheless inter- 
esting to compare the quadratic part of t^ M with 
and i M „. The strongest resemblance occurs when we em- 
ploy our gauge-fixing procedure (that is, we insist that 
hp, v be transverse-traceless) and examine the (0, 0) com- 
ponents of the tensors: 



Too — ^oo — TT^ 1 

.ADM 



i j hij 



dkhijdkhij) 



Q(h 3 ). 



(C6) 



Remarkably, we find that our gauge-fixed Too and too are 
in fact equal to the ADM "Hamiltonian density" tQ UM 
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when working to quadratic order. This is a rather sur- 
prising correspondence, particularly when one considers 
how little our framework has in common with the canoni- 
cal 3+1 approach from which t^ DM arose - Note, however, 
that this equality does not extend to the other compo- 
nents of the tensors: 



KT i0 = -hjkdihjk 



(C7a) 



Kt 



H) 



4 (hjkd 



ihjk — hjkdjhik + hjkdjhki ) (C7b) 



Kt 



A o° M = \ ( hjkdih jk - 2hjkdjh ik ) + 0(h 3 ). (C7c) 



Hence, t^ v , t^„, and do not give rise to equivalent 

local descriptions of gravitational energy-momentum, 
even in transverse-traceless gauge; moreover, even 
though t^ M succeeds in defining a positive gravitational 
energy -density t^ M > as seen by observers at rest 
with respect to the TT-frame, it does not display the 
full (Lorentz-invariant) positivity properties of t^ v : the 
energy-density v^v t/ t^ >M (as seen by an observer moving 
with 4-velocity v^) may be negative, and the energy-flux 
v^t^ M may be spacelike. 

The formulae (C7) also serve as another starting point 
from which to verify that all three tensors define the same 
total momentum CP;: those terms in i i0 an d *io DM which 



do not appear in t,o can be integrated by parts (discard- 
ing a quadratic surface term, as usual) and then van- 
ish due to the gauge condition dihij = 0. Of course, 
these terms do contribute to the total angular momen- 
tum: they correspond to the divergence d a (s^ v a + s„„ Q — 
s a flu )/2 that packages intrinsic spin into the Belinfante 
tensor (22). Accordingly, when one neglects quadratic 
surface terms, one finds that 



'2x\. 



Ad 3 x 



°d 3 , 



2xut 



ADM 0j3 



] U d^ + 0(h d ), 



(C8) 



confirming that t^ v and s a ^ also give the same global 
description of angular momentum as ADM at second or- 
der. 

In summary, whenever general relativity can be ap- 
proximated to quadratic order, and quadratic surface 
terms can be neglected at spatial infinity, our gravita- 
tional energy-momentum tensor t^ v and spin tensor s a ^ u 
provide the same global description of energy, momen- 
tum, and angular momentum as ADM, but localise these 
quantities in a physically sensible fashion, displaying pos- 
itive gravitational energy-density, causal energy-flux, and 
spatial spin. 
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